Analysis of the Knight shift data on Li and Zn substituted YBa2Cu3064 
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The Knight shift data on Li and Zn substituted YBa2Cu306+i! are analysed using an itinerant model 
with short-range antiferromagnetic correlations. The model parameters, which are determined by 
fitting the experimental data on the transverse nuclear relaxation rate T^^ of pure YBa2Cu306+i;, 
are used to calculate the Knight shifts for various nuclei around a nonmagnetic impurity located 
in the Cu02 planes. The calculations are carried out for Li and Zn impurities substituted into 



optimally doped and underdoped YBa2Cu306+a: 
Knight shift measurements on these materials. 



The results are compared with the Li and Y 



PACS Numbers: 74.62. Dh, 76.60.Cq, 71.10.Fd, 74.72.Bk 



I. INTRODUCTION 

The substitution of nonmagnetic impurities into 
cuprates serves as a probe of the electronic correlations in 
these materials. In particular, it has been found that the 
nonmagnetic impurities significantly influence the mag- 
netic correlations in the normal state |lHlo|. The ^^Y 
NMR experiments [0^] find that when Zn is substituted 
into YBa2Cu306+a:, the Knight shifts of the Y sites near 
the impurity are strongly enhanced as the temperature 
is lowered. It has been also found that the Li impurities 
induce changes in the local magnetic environment similar 
to that of Zn . A Curie-like T dependence is found in 
the SQUID measurements of the uniform susceptibility 
for Zn substituted YBa2Cu306-i-a; [§| as well. It has been 
also shown that the Zn impurities induce a Curie-like 
T dependence in the broadening of the ®'^Cu(2) spectral 
line in YBa2Cu306.7 0. The effects of nonmagnetic im- 
purities in YBa2Cu306-i-a; have been also probed by the 
measurements of the ^^Y and ''Li NMR relaxation rates 
Tj~^ and by the inelastic neutron scattering exper- 

iments |ll| Q . The Tj~^ measurements have been also 
carried out for Cu(2) in Zn substituted YBa2Cu307 ||] 
and YBa2Cu408 ^ systems. The effects of the nonmag- 
netic impurities on the magnetic correlations were also 
studied in Lai.85Sro.i5Cu04 with Al impurities pO[ |. 

Numerous theoretical studies p"5|-|2l| have been carried 
out for exploring the effects of the nonmagnetic impu- 
rities on the magnetic properties of the cuprates. The 
effects of a nonmagnetic impurity embedded into the 2D 
Heisenberg lattice p^Jl^ as well as in gapped Heisen- 
berg antiferromagnets |lq| have been studied. The exact 
diagonalization calculations have been used to study the 
effects of a nonmagnetic impurity in the t-J model 
Theoretical studies have been also carried out for the un- 
derdoped YBa2Cu3 06-i-a;, where there is a spin pseudo- 
gap . Calculations of the neutron scattering intensity 
for Zn substituted YBa2Cu306+2; have been also carried 
out |9|||]. 

In this paper, a simple framework pO| of one non- 
magnetic impurity embedded into the lattice of the two- 



dimensional (2D) Hubbard model will be used to make 
comparisons with the ''Li and ®®Y Knight shift measure- 
ments on the Li and Zn substituted YBa2Cu306+a;. The 
purpose of this paper is to study the effects of the non- 
magnetic impurities on the magnetic correlations in the 
normal state of the cuprates. In the following, the 2D 
Hubbard model and the random-phase approximation 
(RPA) will be used for modelling the antiferromagnetic 
correlations in pure YBa2Cu306+2;. The 2D Hubbard 
model is defined by 

H = -t ('^L^J'^ + S^^'^) + ^ 5Z^iT^*TcIiCU 

-l^Y4a'^^'^- (1) 
i.cr 

Here Cia [c]^) annihilates (creates) an electron with spin 
a at site i, t is the near-neighbor hopping matrix element, 
U is the onsite Coulomb repulsion, and /i is the chemical 
potential. The model parameters for the pure system will 
be determined by fitting the experimental data p^-|2^] 



on the transverse nuclear relaxation rate T2 of pure 

YBa2Cu306+a;- 

Next, a static extended impurity potential will be used 
to incorporate the effects of one nonmagnetic impurity 
embedded into the 2D Hubbard model. There is con- 
siderable support for using an extended impurity poten- 
tial for modelling the effective impurity potential . 
Here, it will be assumed that the effective interaction be- 
tween an impurity located at site Tq and a quasiparticle 
at r can be approximated by the static potential 



Kff(ro, r) - VoS{ro, r) + V <5(r, ro + Pa) 



4 
ct=l 



(2) 



which has a range of one lattice spacing. In Eq. (g), a 
sums over the four near-neighbor sites of the impurity as 
illustrated in Fig. 1. Using this simple form of Vcs, the 
Knight shift of various nuclei around the impurity will 
be calculated, and the results will be compared with the 
experimental data. 
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In the following, since the Zn and Li impurities are con- 
sidered to have closed outer electronic shells with S = 0, 
the onsite component Vq will be set to a large negative 
value. On the other hand, the near-neighbor component 
Vi will be used as a free model parameter. Clearly, this 
form of Vcff, Eq- (§)7 is a simple approximation. For in- 
stance, the range of VeS could be more than one lattice 
spacing. Furthermore, Vcs could have scattering in the 
magnetic channel In fact, any magnetic scattering 
could drastically change the results presented here. Nev- 
ertheless, it is interesting to explore the consequences of 
this form of VcS on the Knight shifts within this simple 
framework. Within this model at the level of RPA, the 
fitting of the ^Li and Knight shift data will require 
that Vi is weakly attractive. 

The comparisons will be first carried out with the 
^Li Knight shift, "^K, measurements in Li substituted 
YBasCuaOg.QT by Bobroff et al. §. Here, K will be 
calculated for various values of Vi. It will be seen that 
the temperature dependence of can be fitted by using 
a Vi « — 0.15t. The effects of the nonmagnetic impurity 
on the and the ^^Cu(2) Knight shifts for the sites 
near the impurity will also be given for this compound. 

According to this simple model , the anomalous T 
dependence observed in the Knight shift data is due to 
the enhancement of the antiferromagnetic correlations in 
the local environment of the impurity. In this model, the 
scattering of the antiferromagnetic spin fluctuations by 
the impurity with large momentum transfers near 2ki?, 
where k^? is the fermi wave vector, plays a particularly 
important role [ pOpl| ] . These scatterings locally enhance 
the antiferromagnetic correlations, and, in addition, al- 
low for the Knight shift of the nuclear sites near the im- 
purity to be coupled to the locally enhanced antiferro- 
magnetic correlations. 

Because of the presence of the spin pseudogap in un- 
derdoped YBa2Cu306+a:, it is indeed difficult to extend 
this analysis to this compound. However, the Knight 
shift experiments on underdoped YBa2Cu306+2; yielded 
interesting results on the real-space structure of the in- 
duced magnetic correlations around the impurity. These 
measurements have been carried out for ^Li in Li sub- 
stituted YBa2Cu306.6 Pi and for ^^Y in Zn or Li sub- 
stituted YBa2Cu3 6.64 For ^^Y, two magnetic res- 
onance satellites in addition to the main line are found 
in the presence of the Zn or Li impurities in underdoped 
YBa2Cu306-(-a;. These satellites have been identified as 
belonging to the first and the second neighbour Y sites of 
the impurity. The availability of the Knight shift data on 
^Li and the two ®^Y satellites gives valuable information 
on the real-space structure of the magnetic correlations 
in the local environment of the impurity |^ . The ''Li and 
^^Y Knight shift measurements by Mahajan et al. [Q^] 
and the data on the Knight shift of ^3Cu(2), '^^Kc, by 
Julien et al. Q point out that a staggered polarisation 
cloud forms around the impurity when a uniform mag- 



netic field is applied. Hence, here, using simple assump- 
tions, the results of the calculations will be compared 
with the ^Li and ®^Y Knight shifts in order to see whether 
this simple model has some of the features observed ex- 
perimentally. Results on ^'^Kc will also be presented. It 
is important to keep in mind that the purpose in this 
paper will not be to present a theory of the magnetic 
correlations in the pseudogap state within the presence 
of nonmagnetic impurities. This is clearly beyond the 
simple model used here. Rather, the purpose is to study 
the real-space structure of the deformations induced by 
the impurity by making simple assumptions. It would 
have been more desirable to carry out such a study for 
optimally doped YBa2Cu307, where this model is more 
applicable, but in that case the ^^Y lines are not resolved. 
Here, the model parameters determined by fitting the 
Tj"^ data on pure YBa2Cu3 6.63 will be used to calcu- 
late the ^Li and *^Y Knight shifts. Initially, the Knight 
shifts will be calculated by neglecting the opening of the 
pseudogap in the q — > component of the magnetic sus- 
ceptibility of the pure system. Even in this case, it will 
be found that for Vi between — 0.125t and —Q.15t, the 
real-space pattern of the magnetic correlations around 
the impurity is similar to what is seen experimentally. 
Later on, a pseudo gap in the q — > component of the 
magnetic susceptibility will be introduced artificially. In 
this case, comparisons with the data require Vi « — O.lSi. 
It has to be re-emphasized that the results on the under- 
doped YBa2Cu306-i-2; needs to be interpreted cautiously. 
The treatment of the pseudogap is clearly not rigorous. 
An especially important point is that the presence of the 
pseudogap can introduce new physics which is quite dif- 
ferent than the simple scenario discussed here using a 
weak-coupling model. 

In Section II below, the model will be introduced. Af- 
ter comparisons with the experimental data in Sections 
HI. A and III.B, the pattern of the magnetic correlations 
around the impurity will be shown in Section III.C. The 
implications of these calculations and the role of the an- 
tiferromagnetic correlations will be discussed in Section 
IV. In addition, here, the dependence of the results on 
the effective bandwidth will be studied, and the hyperfine 
couplings used in the calculations will be compared with 
the experimental estimates. In Section V, the summary 
and the conclusions will be given. 

II. MODEL 

Here, the framework for calculating the Knight shifts 
will be introduced. In Section II. A, the modelling of the 
magnetic susceptibility for pure YBa2Cu306+2: will be 
presented. The model parameters, which will be used to 
parameterise the antiferromagnetic correlations of pure 
YBa2Cu306+a;7 wiU be determined by fitting the experi- 
mental data on the transverse nuclear relaxation rate 
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of YBasCuaOe+a; [g2|-g5|. In Section II.B, the static ex- 
tended impurity potential used for approximating the in- 
teraction between the impurity and the electrons will be 
discussed. In the following section, the calculation of the 
magnetic susceptibility for the case of one impurity will 
be presented. In Section II. D, the hyperfine interactions 
for the ^Li and ^^Y nuclear spins will be introduced. 



A. Parameterisation of x(q, for pure YBa2Cu306+a; 

In this paper, it will be assumed that the magnetic 
correlations of pure YBa2Cu306+2; can be approximated 
by the two-dimensional Hubbard model. The magnetic 
susceptibility of the pure system is defined by 



where 



Xpure (q) / dT{m (q,T)m+(q,0)). 



(3) 



Here, r is the Matsubara time, m+(q, 0) ~ 
^^^^^EpCp+qTCpi, m-{q,T) = e^^m,- {q,0)e^"'' , and 
m-(q,0) = A^^^/^EpCp+qiCpt. The RPA form for 

\pure IS 



Xpure (q) 



(4) 



Here, U is the renormalized Coulomb repulsion and the 
static Lindhard susceptibility Xq (q) is given by 



N 1 /(£p+q) ^ /(gp) 

^ ^p ~ ^p+q 



(5) 



where /(Cp) is the fermi factor and Sp — —2t{cospx + 
cospy) ~ fi. Using this form for Xpureio) with the renor- 
malised values of U and the hopping matrix element t, 
the longitudinal and the transverse relaxation rates T^^ 
and T^'^ for YBaaCusOr have been calculated |9|,^. 
In this model, the renormalization of U is due to the 
particle-particle correlations and the single-particle self- 
energy corrections, and the renormalization of the band- 
width is due to the Coulomb correlations |^ . The value 
of U will be determined by fitting the experimental data 
on the ^'^Cu(2) transverse relaxation rate in pure 
YBa2Cu306-i-a;. In the following, the effective bandwidth 
W will be taken to be 1 eV as in Refs. |29 3^. Later, in 
Section IV. C the dependence of the results on W will be 
studied. 

The transverse nuclear relation rate of ^'^Cu(2) for 
an orienting field along the c-axis is given by l22| 



0.69 r 1 



32r 



{^Ei^(q)i'4-e(q) 



^^\M^)\^Xpure{cL] 



(6) 



A(q) = Ac + 2B{cos qx + cos qy) 



(7) 



Here, the following hyperfine interaction has been 
assumed for the ^'^Cu nuclear spin at site i. 



(8) 



5=1 



where 5 sums over the four nearest neighbours of site i. 
For the hyperfine couphngs, Ac — — 4i? = —2.45 x 10^^^ 
erg, corresponding to 164 kOe//iB, will be used as in Ref. 
|||. In fitting the T dependence of T'^ with the RPA 
form of Eq. (^), U will be adjusted as T is varied. Clearly, 
this is a simple, approximate procedure for determining 
the strength of the antiferromagnetic correlations in the 
pure system. In fact, the model is so simple that elec- 
tron filling {n) will be taken to be 0.86 for both the un- 
derdoped and optimally doped systems. Hence, the only 
feature which differentiates between the underdoped and 
the optimally doped systems in this model is the strength 
of the antiferromagnetic correlations of the pure system 
which is set by adjusting U to fit the Tj^^ data. As dis- 
cussed in Rcf. the antiferromagnetic correlations of 
the impure system can get significantly enhanced by the 
impurity scattering as compared to those in the pure sys- 
tem. 

The fitting of Tj"^ will be first carried out for opti- 
mally doped YBa2Cu306+2;. The transverse relaxation 
rate T^^ of ^^Cu(2) has been measured by Pennington 
and Shchter [H in YBasCuaOr, and by Imai et al. H] 
in YBa2Cu306.9. The data by Imai et al. are repre- 
sented by the open circles in Fig. 2(a). Here, it is seen 
that increases linearly as T is lowered from 300 K to 
100 K. The dashed line in this figure has been obtained 
by using the RPA form of Eq. (^) in Eq. (||) and by ad- 
justing U to fit the experimental data. The resultant 
values of U{T) are given by the dashed line in Fig. 2(b). 
Here one observes that as T is lowered from 300 K to 
100 K, U{T) decreases by 5%. Note also that the dashed 
curve in Fig. 2(a) has been extrapolated up to 400 K, 
since in Section III the Knight shifts will be calculated 
up to this temperature. 

The filled circles in Fig. 2(b) represent the T^^"^ data 
for ^'^Cu(2) in YBa2Cu3 06.63 measured by Takigawa [ p5| . 
Here it is seen that Tj""'^ varies almost linearly for T be- 
tween 160 K and 300 K. Below 160 K, T^^ saturates. 
Because of this, the fitting of T^^ for T < 160 K will 
require more attention. This saturation might be due to 
two reasons: (1) the suppression of Xpure{<\ ~ 0) because 
of the opening of the magnetic pseudogap, or (2) the 
saturation of Xpure (q) for q ~ [tt^-k). In the first case, 
XpMre(q ~ (tJ", Ti")) would contiuuc to increase as T is low- 
ered and the saturation of Tj"^ would be due directly 
to the suppression of Xp«re(q ~ 0). In the other case, 
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the saturation of Tj" would be due to the saturation of 
the antiferromagnetic correlations. In order to take into 
account these two possibilities, the fitting of Tj"^ below 
160 K will be done in two ways. In the first case, the 
saturation of T^^ will be attributed to the suppression 
of Xpureio. ^ 0) and it will be assumed that without this 
effect Tj"^ would have continued to increase linearly as 
T is lowered below 160 K. Hence, in this case, Tj"^ will 
be extrapolated for T < 160 K as shown by the solid line 
in Fig. 2(a). The resultant values of U{T) are shown by 
the solid curve in Fig. 2(b). In the second case, the Tj"^ 
data below 160 K will be fitted directly as shown by the 
dotted curve in Fig. 2(a). The resultant U{T) is shown 
by the dotted curve in Fig. 2(b). By comparing the solid 
and the dotted curves for T < 160 K in Figs. 2(a) and 
(b) , it is seen that small changes in U produce significant 
changes in Tj"^. This is because the system is close to a 
magnetic instability with a large Stoner enhancement for 
q - (7r,7r). 

At this point, it is useful to see the strength of the an- 
tiferromagnetic correlations required for fitting the Tj"^ 
data. The dashed curve in Fig. 3(a) shows Xp«re(q) 
versus q obtained by fitting the Tj"^ data on pure 
YBa2Cu306.9 at 100 K. The solid and the dotted curves 
represent Xpure (q) obtained with the first and the second 
scenarios, respectively, for fitting Tj"^ of YBa2Cu3 06.63 
at 100 K For the form of Xp«re(q) used here, Eqs. (^) 
and (^), the peak in Xpure(q) occurs at an incommensu- 
rate wave vector Q* away from Q — (tt, tt). The neutron 
scattering experiments [ |33| , p^ on YBa2Cu3 06+2; find 
that in the superconducting state, the spin-fluctuation 
spectral weight Imxp„re(q, i^) peaks away from (7r,7r), 
but with an incommensuration which is more than what 
is observed in Fig. 3(a). The form of Xpure(q) used here 
is clearly a crude approximation. The results shown be- 
low will not depend sensitively on the amount of the in- 
commensuration, rather they will depend on the total 
integrated weight in the q ~ (tt, tt) region of the q-space. 
Further results on the nature of the antiferromagnetic 
correlations are shown in Fig. 3(b). The dashed line in 
Fig. 3(b) shows Xptire(Q*) versus T for optimally doped 
YBa2Cu306+a;- Also shown in this figure by the solid 
line are the results for x(Q*) of f^ic impure system with 
0.5% randomly distributed impurities which will be dis- 
cussed below in Section IV. A. Fig. 3(c) shows similar 
results for x(Q*) of underdoped YBa2Cu3 06-i-x- Here, 
it is seen that x(Q*) g^ts significantly enhanced by the 
substitution of the nonmagnetic impurities, especially for 
underdoped YBa2Cu3 06+a:. By comparing Fig. 2(a) with 
Figs. 3(b) and (c), it is also seen that the temperature 
dependence of Tj"^ and Xpure{Q*) are closely related, as 
expected. 



B. Effective Impurity Potential 

The effects of one nonmagnetic impurity located at site 
Tq will be taken into account by adding to Eq. (|l|) the 
term 

^ Voff(ro,r,)ctc,o. (9) 

where VJ;ff(ro,r) is given by Eq. (||). The importance of 
using an extended impurity potential has been previously 
noted [^,1^. In Ref. it has been pointed out that 
the Coulomb correlations of the host could cause the ex- 
tended nature of the effective interaction. The fact that 
the defects generated by the electron irradiation of the 
samples have similar effects on the Tc suppression ]3^ ] 
and on the magnetic properties js^ also supports using 
an impurity potential with only a potential scattering 
term. 

In the following, the onsite component of the impu- 
rity potential, Vb, will be set to a large negative value, 

— loot, in order to model the closed electronic shell of the 
nonmagnetic impurity. As long as |Vo| has a large value, 
its exact magnitude or its sign does not play an impor- 
tant role. For instance, using Vq = — 200t instead of 

— loot does not change the results. The near-neighbour 
component Vi, on the other hand, will be used as a free 
parameter in fitting the Knight shift measurements on Li 
and Zn substituted YBa2Cu3 06-(-x- In the next section, 
the effects of this static extended impurity potential on 
the magnetic susceptibility will be calculated. 

C. Magnetic susceptibility within the presence of 
one nonmagnetic impurity 

The calculation of x within the presence of a non- 
magnetic impurity used here follows that given in 
Ref. First, the effects of the impurity on the single- 
particle Green's function will be calculated. The single- 
particle Green's function in Matsubara frequency space 
is defined by 

G(r„r,,ic^„) = -^''dre*'^"^(c,,(T)c],(0)), (10) 

where tOn — (2n + l)7rr. For the pure system with [/ = 0, 
one has 

Go(r„r„io;„) = ^ ^ e'P-('^--"-^)Go(p, zc^„) (11) 
p 

where 

Go{p,iu;n) = -. ^- . (12) 

If an impurity is introduced at site Tq, then one gets 
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G{r,r',iujn) = Go(r, r', za;„) 

+ ^ Go(r,r",zc^„)T(r",r"',zc^„)G(r"',r',ic^„), (13) 

r" .r'" 

where the T-matrix for the impurity scattering is given 

by 

T(r,r',zw„) = ,5(r,r')ycff(ro,r) 

+ J2 ^off(ro, r)Go(r, r", ia;„)T(r", r', za;„). (14) 

r" 

The calculation of G(r, r', iujn) is illustrated diagrammat- 
ically in Fig. 4(a). 

When the translational invariance is broken, the mag- 
netic susceptibility is defined in real space as 



Jo 



(TO-(r,r)m+(r',0)), (15) 



where ■m'^{r) — cl^Cri, and TO^(r) — cl^Cri- The effects 
of a single impurity will be first calculated for U = 0, giv- 
ing the irreducible susceptibility xo(r, r', iwm), and then 
the effects of the Coulomb correlations will be included. 
The diagrams representing the effects of the impurity on 
Xo are shown in Figs. 4(b) and (c). Both the self-energy 
and the vertex corrections need to be included |^^, and 
the resulting expression for xo(r, r', it^m) is given by 



G(r, r',iw„ + icj^) 



xGo(r',r, iw„) + Go(r, r',iw„ + iuj,n)G{r' ,r,iujn) 
-T^ ^ Go(r, ri,ia;„ -t-iLj„i) 

iuJn I"i,r2,r3,r4 

xGo(r2,r',iu;„ -I- ia;,„)T(ri, r2, iw„ + iujm) 
xT(r3, r4, iu;„)Go(r3, r, zw„)Go(r', r4, iw„). (16) 

Note that the irreducible impurity-scattering vertex has 
been used in calculating the vertex corrections to xo in- 
stead of the reducible one. This is necessary in order to 
prevent double counting, since here the effects of only 
one impurity are calculated. 

Next, the Coulomb correlations are included by solving 
the RPA equation 

X(r, r', iujm)^ Xo(i", r', iuj^) 

+U J2 Xoir, r", zc^„)x(r", r', iuj„,) (17) 

r" 

for x(r, r', iujm)- Here x is calculated for an A^-site square 
lattice with one nonmagnetic impurity located at the cen- 
ter using periodic boundary conditions. The calculations 
are carried out on sufficiently large lattices so that the 
finite size effects are small. In Ref. the finite size 
effects on x have been studied in detail. 



The Knight shift is determined by xI^jT') = 
x(r,r',*Wm = 0). In the following section, the hyperfine 
interactions, which relate x(r, r') to the Knight shifts, 
will be discussed for ^Li and *^Y. 



D. Hyperfine interactions for ^Li and 

The hyperfine interaction between the ^Li nuclear spin 
with ''1 = 3/2 substituted into a Cu(2) site and the elec- 
tronic spins will be modelled by the following hyperfine 
coupling 



(18) 



5=1 



Here S sums over the four Cu(2) sites neighbouring the 
Li impurity. The nature of this interaction is similar to 
that of the transferred hyperfine coupling of the ^^Cu(2) 
nuclear spins to the electronic spins at the neighbouring 
Cu(2) sites given by the second term in Eq. (^) |Q. In 
Section III, the magnitude of the ''Li hyperfine coupling 
G will be taken to be 1.8 x 10~^° erg, corresponding to 
0.85 kOe/nB- This choice for the value of G will be dis- 
cussed later in Section IV. C. As a result of this coupling, 
Eq. (|l8|) , the Knight shift of ''Li is (see Appendix I) 



G4fc(r=(l,0)), 



where fc(r) is defined by 



fc(r)=^x(r,r'), 



(19) 



(20) 



and 7e and "^^n are the gyromagnetic ratios of the electron 
and the ^Li nuclear spin. The factor of four in Eq. (|l|) 
is because the Li impurity has four Cu(2) neighbours. 

The hyperfine interaction for the ^^Y nuclear spin with 
89/= 1/2 is 



''I -So 



(21) 



where a sums over the eight Cu(2) sites neighbouring the 
^^Y nuclear spin. In Section III, the hyperfine interaction 
D will be assumed to be —2.2 x 10~^° erg, corresponding 
to —1.0 kOe/nB, in order to fit the normal state value of 
^^K in pure YBa2Cu306+x using a bandwidth of = 1 
eV. Later in Section IV. C, the dependence of the results 
on W and on the values of the hyperfine couplings G and 
D will be discussed. Here, it will also be assumed that 
the ^^Y hyperfine couplings do not change upon the sub- 
stitution of the nonmagnetic impurity. Then the Knight 
shift for at a particular Y lattice site is given by 
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89-, ,D 2^ '(■■»). 

In ; 



(22) 



where Tq, are the locations of the eight nearest-neighbour 
Cu(2) sites of the Y site in the Cu02 bilayer. For dilute 
Li or Zn impurities, the ^^Y nuclear spin will be affected 
by an impurity located in one of the layers of the Cu02 
bilayer in YBa2Cu30g+a;. Here, it will be assumed that 
the changes induced in ^'^K are due only to the changes in 
the magnetic correlations in the layer which contains the 
impurity as it was done in Ref. [Q. Also, in the following, 
the quantity of interest for *^Y will be the change in the 
Knight shift of '^^Y due to the presence of the impurity. 
Hence, the calculations will be carried out for 



2 



89 



In 



D 



A:(r„) - 4x 



pure 



(23) 



where a sums over the four Cu(2) sites neighbouring 
the impurity in the layer which contains the impurity. 
In Eq. (H), 

Xpure is the uniform static susceptibility of 
the pure system. The value of ^^K will depend strongly 
on the location of the ^^Y nuclear spin with respect to 
the impurity. In Fig. 5, a sketch of the various Y sites 
with respect to the impurity site is given. Note that the 
measurement of ^^K at the Y sites near the impurity, 
along with that of K, determines the real-space struc- 
ture of fc(r) in the environment of the impurity. For in- 
stance, ^'^K for Y(l) is determined by fc(r = (1,0)) and 
k{v = (1,1)), while '^K is set by fc(r = (1,0)). In Sec- 
tion HI, the results of the calculations on K and ^'^K 
for various Y sites will be shown and compared with the 
experimental data. 



III. RESULTS OF THE CALCULATIONS 

The comparison of the results with the experimental 
data will be carried out first for the optimally doped and 
then for the underdoped YBa2Cu306+a:. 

A. Results on optimally doped YBa2Cu306+a; 

The filled circles in Fig. 6(a) denote the experimental 
data by Bobroff et al. on the ^Li Knight shift K for 
YBa2Cu306.97 with dilute Li impurities. The curves in 
this figure represent the results of the calculations of K 
for various values of Vi as indicated next to the curves. 
The model parameters other than Vi were already set by 
fitting the Tj"^ data on pure YBa2Cu30g.g by Imai et al. 

as discussed in Section II. A. Here, it is seen that for 
Vi ~ — 0.15i the calculated values of K are in agreement 
with the experimental data. It is useful to compare these 
results with what is expected for the ^Li Knight shift 



if the magnetic correlations had not changed upon the 
substitution of the impurity and had remained the same 
as in the pure material. In Fig. 6(b), Kq defined by 
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In 



(24) 



where Xpure is the uniform susceptibility of the pure sys- 
tem, is plotted as a function of the temperature. The 
solid line has been obtained from Xp«re(q — * 0) given by 
Eq. (Q) for optimally doped YBa2Cu306+i. The dashed 
line in Fig. 6(b) has been obtained by rescaling the ^^Y 
Knight shift ss/iT of pure YBa2Cu307 ISSl with 
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The factor of 1/2 in this expression is because Y has eight 
Cu(2) near neighbours while Li has four. In these figures, 
it is also seen that for Vi = 0, corresponding to an onsite 
impurity potential in the unitary limit, K is enhanced 
with respect to ^A'o but this enhancement is not sufficient 
to explain the experimental data. If positive values of Vi 
are used, then K gets suppressed with respect to the 
Vi = case. 

The stoichiometry of the samples is a factor which 
could affect the comparisons with the data. Note that 
U{T) used in calculating K for YBa2Cu306. 97 was 
determined by fitting the T2~^ data on YBa2Cu306.9. 
For YBa2Cu307, T^"^ was measured at 100 K, and its 
value is 7.7 ± 0.6 msec~^ compared to ^ 9.5 msec~^ for 
YBa2Cu306.9. The value of 7.7 msec^^ for T2"~^ requires 
a slightly smaller value for U and, consequently, the fit- 
ting of ^if is done using larger values of Vi. If one as- 
sumes that the T dependence of is given by a line 
which passes through 7.7 msec~^ at 100 K and which is 
parallel to the dashed line in Fig. 2(a), then the K data 
can be fitted by using Vi ~ — 0.275t. 

Based on these results and considering the simplicity of 
the model, the conclusions which can be reached are lim- 
ited. Probably, the best thing to say is that the analysis 
of the ^ K data on optimally doped YBa2Cu306+2; within 
this model at the level of RPA requires an impurity po- 
tential which is weakly attractive at the nearest-neighbor 
Cu(2) sites. 

The results of the calculations for ^'^K are shown in 
Fig. 7(a). Here, the temperature dependence of is 
shown for the first four nearest-neighbor Y sites which 
are indicated in Fig. 5. In obtaining these results, 
Vi ~ — 0.15i was used. Note that in Fig. 7(a) the ex- 
pected values of the Knight shifts are close to each other. 
Experimentally, for optimally doped YBa2Cu306+2;, the 
expected resonance lines for the ^^Y nuclei near the Zn 
impurity are not resolved. 

Next, results on how the nonmagnetic impurity affects 
the ^■^Cu(2) Knight shifts in this model are presented. 
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The ^'^Cu(2) Knight shift at site with the orienting 
magnetic field H||c is given by 
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where a sums over the four Cu(2) neighbors of r^. Figure 
7(b) shows the T dependence oi^^K^ivi) for Vi = -0.15t 
for the first four neighboring Cu(2) sites of the impurity, 
which correspond to = (1,0), (1,1), (2,0), and (2,1). 
For the pure system, ^^K^ vanishes since Ac + 4i? = 0. 
Here, it is seen that ^^Kc has a staggered pattern; it is 
negative at the sublattice of the (1,0) site and it is posi- 
tive at the sublattice of the impurity. This spatial pattern 
oi^^Kc{vi) is due to the staggered pattern of A;(r), which 
will be discussed in Sections III.C and IV. A. The mea- 
surements oi^^Kc have been carried out in YBa2Cu306.7 
by Juhen et al. 0. While the individual ^^Cu(2) lines 
are not resolved, it has been observed that the broaden- 
ing of the linewidth has a Curie-like T dependence. This 
has been attributed to the development of a straggered 
polarization cloud around the impurity j^]. 

B. Results on underdoped YBa2Cu306+a: 

The analysis of the Knight shift data on the under- 
doped YBa2Cu306+x is considerably more complicated 
because of the magnetic pseudogap. This is so because 
the origin of the magnetic pseudogap in this material 
is currently an unresolved issue. Hence, it is necessary 
to note that the purpose here is not to develop a the- 
ory for the magnetic susceptibility of the underdoped 
YBa2Cu30g+a; with impurities. Instead, the purpose 
is to explore any possible role of the antiferromagnetic 
correlations in producing the anomalous T dependence 
of the Knight shift data in underdoped YBa2Cu306+i 
within the presence of nonmagnetic impurities. The em- 
phasis here will be on the real-space structure of the mag- 
netic correlations around the impurity. 

The ^Li Knight shift measurements by Bobroff et al. ||] 
on underdoped YBa2Cu306+a; with Li impurities found 
that K increases rapidly as T is lowered. In Fig. 8(a), 
the experimental data on 5'' K defined by 



K Kq 



(27) 



are compared with the results of the calculations. The 
comparisons are carried out for K rather than for K 
in order to compensate for the effect of the magnetic 
pseudogap which is not taken into account here. The 
data points shown by the filled circles in Fig. 8(a) were 
obtained by subtracting from the '^K data of Bobroff et 
al. M the quantity 
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if (main). 



(28) 



where if (main) is the Knight shift of the main ^^Y 
spectral line measured by Mahajan et al. [Q in Zn sub- 
stituted YBa2Cu306.6- The solid points in Fig. 8(b) 
indicate Kq estimated this way. The dashed curve in 
Fig. 8(b) shows the estimate for Kq obtained by using 
the measurements on pure YBa2Cu3 06.63 [§8|. In 
the Zn substituted samples, *^i4r(main) is slightly shifted 
with respect to ^^K of the pure samples. This might be 
due to the difference in the stoichiometry of the Zn sub- 
stituted and the pure samples, as discussed in Ref. Q. 

The results of the calculations on S'' K are represented 
by the various curves in Fig. 8(a). The term subtracted 
from K to obtain these results is shown by the solid line 
in Fig. 8(b). Note that the sohd line in Fig. 8(b) lies 
slightly above the solid line in Fig. 6(b), since larger val- 
ues of U were used for the underdoped YBa2Cu3 06+x. 
These calculations are shown for four different values 
of Vi: -Q.l2bt (short-dashed), -0.15t (solid), -0.175i 
(long-dashed), and 0.0 (dot-dashed). These results have 
been obtained by using U{T) represented by the solid 
curve in Fig. 2(b). Hence, they correspond to the first 
scenario described above where Tj"^ has been linearly 
extrapolated for T < 160 K. For each value of Vi, the 
calculations of S'^K have been repeated according to the 
second scenario, where the saturation of T^^ below 160 
K is attributed to the saturation of the antiferromagnetic 
correlations. In this case, the values of U{T) shown by 
the dotted curve in Fig. 2(b) were used. In Fig. 8(a), 
these results are shown by the dotted curves for each 
value of Vi as indicated. Comparing the results of the 
calculations with the experimental data, one sees that 
the values of Vi between — 0.125i and —0.15t would pro- 
duce a fit of the experimental data on J^iC. 

Next, the results for *^Y are discussed. The experi- 
ments by Mahajan et al. found that Zn substitution 
in underdoped YBa2Cu306+2: strongly modifies the nu- 
clear resonance spectrum of ^^Y. In this case, in addition 
to the main ^^Y resonance line, two satellite peaks are ob- 
served. These outer and middle resonance lines have been 
identified as belonging to the ^^Y nuclear spins which are 
respectively the first and the second near-neighbours of 
the Zn impurity, Y(l) and Y(2) in Fig. 5. The fact that 
the magnetic resonance spectrum of ^^Y in Li substituted 
YBa2Cu306+a; is nearly the same as in Zn substituted 
samples confirms this identification ||^. In the following, 
the data on ^^K will be compared with the results of the 
calculations which were outlined in Section III. However, 
the comparisons will not be carried out directly for ^^K. 
Rather, the shift of the satellite lines with respect to the 
main line will be used in order to compensate for the 
opening of the pseudogap in the pure system. Hence, the 
comparisons will be carried out for the following quanti- 
ties, 



^^(outer) X(outer) X(main) 
^if(middle) ii:(middle) -^^ K{mam). 



(29) 
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Since *^ii'(main) in Zn substituted material nearly fol- 
lows the T dependence of ^^K in the pure system, 
5^^K{outei) and (5^^ if (middle) represent the change in 
the Knight shifts of the nuclear spins at the Y(l) and 
Y(2) sites, respectively, due to the substitution of the Zn 
impurity. 

In Figures 9(a) and (b), the points represent the ex- 
perimental data on —S^^K{outeT) and —i5^^ if (middle), 
respectively, for underdoped YBa2Cu306+a; measured by 
Mahajan et al. 0,^- In these figures, the data points are 
shown for two different orientations of the magnetic field 
H. It is estimated that the ^^Y hyperfine coupling has an 
anisotropy of about 15%, which could be the cause of the 
anisotropy of the data seen in Figs. 9(a) and (b) fl^. 
In Fig. 9(a), the error bars shown for the lowest and the 
highest temperature measurements are representative of 
the error bars for the other data points. In Fig. 9(b), 
the data points arc shown up to 130 K, since at higher 
temperatures the middle satellite is not resolved. 

The curves in Figs. 9(a) and (b) represent the results 
of the calculations for —J^^if (outer) and — if (middle), 
respectively. These calculations have been carried out for 
the same values of Vi as in Fig. 8(a) for S'^ K [Q. For 
T < 160 K, these calculations were carried out in two 
ways by using U{T) shown by the solid and the dotted 
curves in Fig. 2(b). In Figs. 9(a) and (b), it is seen that 
for Vi between — 0.125t and — O.lSt, the data can be fit- 
ted. Finally, in Fig. 10 the T dependence of ^^ifc(i"i) 
in underdoped YBa2Cu306+a; is shown for various Cu(2) 
sites near the impurity. Here, Vi = —0.15t was used. In 
order to interpret these, in the next section the real-space 
structure of the changes induced in the magnetic corre- 
lations by the impurity will be shown. The ^^K^ mea- 
surements in Zn substituted YBa2Cu306.7 imply that the 
polarisation around the impurity is staggered [[7| . 

C. Pattern of the magnetic correlations around the 
impurity 

In Fig. 11(a), k{r) defined by Eq. (|^), is shown as 
a function of the lattice distance r = |r| away from the 
impurity at T = 100 K. Here, r is given in units of the lat- 
tice constant, and in obtaining these results Vi = —0.15t 
was used. The horizontal long-dashed line denotes what 
is expected for k{r) within this model for pure opti- 
mally doped YBa2Cu3 06+2;, and the open circles repre- 
sent fc(r) when an impurity is introduced into this mate- 
rial at ro = 0. The results for underdoped YBa2Cu306+2; 
are also shown in Fig. 11(a). The dotted horizontal line 
denotes k{r) for pure underdoped YBa2Cu306+2;, while 
the filled circles show fc(r) within the presence of the 
impurity for this material. The results on underdoped 
YBa2Cu306+a; were obtained using the value of U given 
by the solid curve in Fig. 2(b). Note that, for under- 
doped YBa2Cu306+£c in this graph, the physically more 



relevant quantity is the change in k{r) induced by the 
impurity, since the presence of the pseuogap in the uni- 
form susceptibility of the pure system is not included in 
the model. 

The real-space structure of fc(r) is important since it 
determines the Knight shift at various sites. For instance, 
fc(r) at r = (1, 0) determines the value of ^if . On the 
other hand, ^^^if (outer), corresponding to the Y(l) site 
in Fig. 5, is given by the change in 

fc(0,0)-|-2fc(l,0) + fc(l,l) (30) 

induced by the impurity. Similarly, ^*^if (middle), cor- 
responding to the Y(2) site, probes the change induced 
in 

fc(l,0) + fc(l,l) + fc(2,0) + fc(2,l) (31) 

by the impurity. Hence, it is seen that by the fitting of 
S'^K in Fig. 8(a), the value of k{l,0) is fixed. Since it 
is known that k(0,0) vanishes when the impurity is in- 
troduced, the fitting of J^^if (outer) in Fig. 9(a) fixes the 
value of fc(l, 1). In Fig. 11(a), it is seen that while fc(l, 0) 
is enhanced strongly by the impurity, fc(l, 1) is strongly 
suppressed and it even becomes negative. The fact that 
both (5^if and J^^if (outer) can be fitted at the same time 
means that the induced changes in fc(r) of the actual sys- 
tem must have a staggered pattern at the sites (1,0) and 
(1,1). Since fc(l, 0) and fc(l, 1) are now fixed, the fact that 
(5^^if (middle) can also be fitted means that the change 
in fc(2,0) + fc(2, 1) is fixed. In Fig. 11(a), it is seen that 
fc(2, 0) gets suppressed by the impurity and fc(2, 1) gets 
enhanced while the total effect on k{2, 0)+k{2, 1) is small. 
This in turn means that the induced changes in k(r) of 
the actual system must also be staggered at sites (2,0) 
and (2,1). Hence, the fits obtained in Figs. 8(a), 9(a) 
and (b) imply that k(r) has a staggered pattern in the 
vicinity of the impurity. This conclusion was reached by 
the analysis of Mahajan et al. [0^]. Here it is seen that 
this simple model produces this pattern as well. 

It is useful to compare the staggered pattern of k{r) 
with the pattern of x(r, r) near the impurity. Fig. 11(b) 
shows x(r,r) versus r plotted in the same way as fc(r) 
is plotted in Fig. 11(a). Here one sees that x(r,r) is 
strongly modified near the impurity. It is also seen that 
x(r, r) does not have a staggered pattern as, for instance, 
one would have had for antiferromagnetically ordered 
spins. The comparisons shown in Figs. 6, 8 and 9 mean 
that the experimental data on ^if and ^^K require a 
staggered pattern for fc(r) in the vicinity of the impurity. 
But a staggered structure for x(r, r) is not necessary. 

In Figs. 11(a) and (b), it is seen that while x(r,r) is 
positive for all r, fc(r) can become negative, for instance 
at r = (1,1)- This is possible because x(r, r') for r = 
(1,0) and r' = (1,1) has a large negative value after 
the impurity is introduced, meaning that the impurity 
is inducing strong antiferromagnetic correlations in its 
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local environment in this model. These points will be 
discussed further in Section IV. A. 

At this point it is also useful to present results on how 
the impurity potential affects the single-particle prop- 
erties. For this purpose, in Fig. 11(c) the electron oc- 
cupation number n{ri) = X]o-('^L'-jo') sites near the 
impurity is shown as a function of the distance r = \ri\ 
from the impurity for various values of Vi/t. Since Vq is 
strongly attractive, the impurity site which is not shown 
in this figure, is nearly doubly occupied, n{ri = 0) « 2.0. 
For Vi = 0, it is seen that n{ri) is enhanced at = (1, 1) 
with respect to its value at — (1,0). Away from the 
impurity, these oscillation in niji) decay, and n{ri) goes 
to its value in the pure system, (n) = 0.86, which is indi- 
cated by the horizontal long-dashed line. Here, it is also 
seen that this structure in n(ri) holds for weak attractive 
values of Vi. But as Vi becomes more attractive, for in- 
stance for Vi — — 0.15i, the pattern in n(rj) changes. In 
this case, n{ri = (1, 0)) gets enhanced over n(ri = (1, 1)). 
Hence, in Fig. 11(c) one observes that the structure in 
n(ri) in the vicinity of the impurity does not necessarily 
reflect the sign of Vi. For instance, a finite attractive 
Vi is required before n{ri — (1,0)) gets enhanced with 
respect to n{ri = (1, 1)). This has to do with the pres- 
ence of the strongly attractive Vq term in the impurity 
potential. 

In Section III.B, it was seen that the data on (J^if 
and S^^K can be fitted reasonably well over the whole 
temperature range by using a simple form for the effec- 
tive impurity potential without any T dependence. On 
the other hand, because of the various approximations 
employed in this model, the quantitative agreement ob- 
tained in the fits might actually be misleading. However, 
what is significant is the fact that S'^K, S^^K{outei) and 
i5^^iC(middle) can be fitted all at the same time, which 
means that the real-space structure of the changes in- 
duced in the magnetic correlations near the impurity ap- 
pears to be described by this simple model. 

IV. DISCUSSION 

A. Role of the antiferromagnetic correlations 

It is useful to discuss the connection between the an- 
tiferromagnetic correlations of the system which is com- 
posed of the impurity and the host, and the anomalous 
enhancement of the Knight shift at sites around the im- 
purity. This connection becomes more clear if the Fourier 
transform xi^^^l') is introduced, 

x(q,q') = E^''''""''''""'^^(^'^')- (32) 

Note that x(r, r') is the susceptibility for an A^-site 
square lattice with one impurity at the center and pe- 
riodic boundary conditions. 



An especially important quantity is the off-diagonal 
susceptibility xi.^^^) = xlq^q' = 0); since its Fourier 
transform with respect to q gives fc(r). Particularly, at 
sites r = (1,0) and (1,1) with respect to the impurity, 
one has 

q 

/c(l,l)=Ecos((7, + g^)x(q,0). (33) 
q 

In Figs. 12(a) and (b), the q dependence of — Xo(q, 0), 
which is for U = 0, and — x(q7 0) ^ire shown at 100 K 
for optimally doped YBa2Cu306+a: obtained by using 
Vi = — 0.15i l^lj]. Here, the (5-function component of 
x(q; q' = 0) at q = has been omitted. In these fig- 
ures, it is seen that — x(qiO) is enhanced with respect to 
— Xo(q, 0), and they both peak at q near (7r,7r). Note 
that for q ^ (tt, tt), the form factors entering Eq. (p^), 
i(cos(7j. + cos qy) ^ —1 and cos{qx + qy) ~ 1. Because 
of the q structure of x(qiO), fc(l,0) gets enhanced and 
fc(l, 1) gets suppressed. Hence, within this model, the 
staggered pattern of /c(r) near the impurity is a conse- 
quence of the peaking of —x{^.^ 0) near (tt, tt). 

The enhancement of — x(qi 0) with respect to 
— Xo(q,0) is understood better if the RPA equation, 
Eq. (|l7|), is written in momentum space, 

x(q, q') = xo(q, q') + E ^o(q' q")x(q", q'), (34) 

q" 

where Xo(q, q') is for U — Q. For the pure system, 
X(q, q' = 0) is given by 

X(q,0) = 7V(5q,oXp«re(q^O) (35) 

where Xpnre is defined by Eq. (^ , and within the presence 
of the impurity x(q, q' = 0) is obtained by solving 

X(q,0) =Xo(q,0) + f/Exo(q",0)x(q",q). (36) 
q" 

This expression shows that when q" ~ (7i',7r), x(q ^ 
(tt, 7r),q' = 0) couples to the antiferromagnetic correla- 
tions determined by x(q, q) with q ^ (tt, tt) of the system 
which is composed of the host and the impurity. This is 
the reason for the strong enhancement of x(q, q' — 0) 
with respect to xo(qi q' = 0). 

It is useful to discuss the physical meaning of the off- 
diagonal susceptibility Xo(q, q') where q 7^ q'. Dur- 
ing the scattering of the spin fluctuations by the im- 
purity potential, the momentum is not conserved and 
Xo(q7q 7^ q') acts as the vertex for the scattering of 
the spin fluctuations by the impurity with Q* = q — q' 
momentum transfers. The peaking of Xo(q5q' — 0) near 
q ~ (tt, tt) means that the scattering of the antiferromag- 
netic spin fluctuations with momentum transfers near 
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(tt, tt) is the dominant scattering process. The effects 
of the scatterings with large momentum transfers on the 
Q = (tt, tt) neutron scattering intensity have been also 
emphasized in Ref. where Xo(q, q') was calculated 
at the lowest order in the impurity potential. Within 
this model, both the Knight shifts and the neutron scat- 
tering experiments point out at the importance of the 
scattering of the particle-hole pairs by the nonmagnetic 
impurity with large momentum transfers near ikp. If 
the ~ 2ki7' scatterings are indeed one of the primary ef- 
fects of the nonmagnetic impurities, then an anomalous 
softening of the phonons at wave vectors ^ 2]ip might 
be observed in YBa2Cu306+2; with dilute Zn impurities 
11- 

The real-space pattern of ^^Kc{ri) seen in Fig. 7(b) 
also reflects the q dependence of x(q, q' = 0). Equa- 
tion ( |2^ ) for ^^Kc{ri) can be rewritten as 

''^c(r.) = i (^) 1 V e^'^ - (A, + 4S7,)x(q, 0), 



(37) 



where 7q = (cosg^; -I- cos(?y)/2. Since x(q7q' = 0) 
most of its weight at q (tt, tt) for which 7q ~ —1, one 

gets 
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In 



N ^ , 

q~(-7r,7r) 



x(q,o). 



(38) 



Here one sees that 4- \B has a large value and 
^'^Kcivi) directly couples to x(q ^ (^, 7r),q' — 0). It has 
been noted in Ref. |^] that the anomalous broadening of 
the ^'^Cu(2) linewidth reflects the staggered polarization 
of the magnetic correlations near the impurity. 

It is important to note that, within this model, the an- 
tiferromagnetic correlations also get enhanced due to the 
presence of the impurity Q . If Ni randomly distributed 
nonmagnetic impurities are considered in the dilute limit, 
then the q-dependent magnetic susceptibility of this sys- 
tem is given by 



x(q,q) - Xp«re(q) 



(39) 



where rii — Ni/N is the impurity concentration. The 
results on x(Q*) for 0.5% nonmagnetic impurities are 
compared with Xpure 

(Q*) in Figs. 3(b) and (c), where it 
is seen that the impure system has stronger antiferromag- 
netic correlations than the pure system. Here, Q* is the 
incommensurate wave vector where x peaks at low tem- 
peratures. The enhancement of x(Q*) over Xpure(Q*) is 
especially significant for the underdoped system. Hence, 
it needs to be noted that the enhancement of the Knight 
shifts near the impurity in this model is a result of the 



couphng to x(q; q) for q ^ {tt, tt) of the whole system 
which is composed of the impurity and the host rather 
than to Xp«re(q ^ (tt, tt)) of the pure host. 

Experimentally, Imx(Q,t^) is the inelastic neu- 
tron scattering spectral weight. The enhancement of 
Imx(Q,w) in the normal state by the impurity scatter- 
ing is clearly seen for dilute Zn impurities in YBa2Cu307 
p2|-p^ . The low frequency part of Im x is also probed 
by the measurement of the longitudinal relaxation rate 
T^^ at sites near the impurity. The measurements of 
T^^ at sites near the impurity have found interesting re- 
sults [0J^,^,01 • The calculations presented here have been 
extended to obtain the ^Li Tj~^ using the Pade approx- 
imation for analytic continuation to the real-frequency 
axis of the results calculated in terms of the Matsubara 
frequencies. While it is difficult to obtain reliable results 
on Imx(Q,w) for general to using the Pade approxima- 
tion, it is possible to obtain control on the calculation 
of T^^, which requires only the lu —> limit. Since the 
analytic continuation procedure requires special atten- 
tion, the calculations of the ^Li T^^ will be presented 
elsewhere 



B. Effects of the pseudogap in underdoped 

YBa2Cu306+. 

The analysis of the data on underdoped YBa2Cu306+a; 
needs to be interpreted carefully because of the presence 
of the pseudogap in this system. The uniform suscep- 
tibility of pure YBa2Cu3 06.63 starts to decrease below 
300 K, with a T dependence proportional to that given 
by the dashed curve in Fig. 8(b). This has been ignored 
in the calculations presented above. In order to explore 
the effects of the pseudogap on the Knight shift results, 
the following simple calculation has been carried out. 
The suppression of the diagonal irreducible susceptibility 
Xo (q, q) for q ^ has been artificially incorporated into 
the model by multiplying xo (q^ q) by the factor 



F(q) = 1 - ae-l'il'/'^' 



(40) 



which has been chosen in order to cause a suppression 
for q ~ 0. After artificially suppressing the diagonal ir- 
reducible susceptibility for q ^ 0, the Fourier transform 
is taken to obtain xo(r,r') of the pure system, which is 
then used in calculating xo(r;r') of the impure system. 
Hence, in this procedure, the diagonal terms Xo(qj q) for 
q ~ of the impure system have the effects of the pseu- 
dogap, but Xo(q, q) for q ~ (7r,7r) and the off-diagonal 
terms Xo(qjq'); where q ^ q', are not affected by this 
artificial opening of the pseudogap. Next, xo(r,r') cal- 
culated with this procedure was used in solving Eq. ( [T^ ) 
for x(r,r') and in determining the Knight shifts. The 
parameter k entering Eq. (^0|) was arbitrarily chosen to 
be tt/2. In addition, a has been chosen such that the 
ratio of 
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1-C/F(0)x^(0) 



(41) 



to Xo (0)/(l~^Xo(0)) is equal to the ratio of^^K in pure 
underdoped YBa2Cu306+x to that in pure optimally 
doped YBa2Cu3 06+a;- For instance, at T = 100 K, this 
condition is satisfied by using a = 0.51. The opening of 
the pseudogap in this way requires slightly larger values 
of Vi for fitting the Knight shift data. In Section III.B, 
the data on S'^K, S^^K{outej:) and ^*^i4r(middle) were 
fitted by using Vi between — 0.125i and — 0.15t. Here, 
if Vi — — 0.15i is used along with the above mentioned 
values of k and a, then at T = 100 K one obtains the fol- 
lowing results for the Knight shifts: S'^K — 650 ppm, 
S^^K{outer) = -180 ppm, and ^^^^^^(middle) = -75 
ppm. On the other hand, at 300 K, it is necessary to use 
a = 0.25 and in this case one obtains: K = 250 ppm, 
(5*^ii'(outer) = —35 ppm, and 5^^i4r(middle) = —20 ppm. 
These values for the Knight shifts are comparable to the 
data seen in Figs. 8(a) and 9. Hence, when the pseudo- 
gap is introduced in Xo(qiq) for q in this artificial 
way, the Knight shift data can be still fitted, but by using 
a slightly larger value for V\. However, it must be kept 
in mind that the way the pseudogap is introduced here 
is not rigorous, and in fact the fitting of the data in this 
way can be considered as supcrfiuous. For this reason, 
the fitting of the data on the underdoped compound will 
not be pursued further. Rather, here it is only pointed 
out that the staggered nature of the induced magnetic 
correlations seen in this model appears to be consistent 
with the measurements of the ^Li and Knight shifts 
in underdoped YBa2Cu3 06+a;- 



C. Dependence of the results on the effective 
bandwidth 

The results seen in Section HI were obtained by us- 
ing an effective bandwith of 1 eV. For a one-band 
model of the cuprates, the bare hopping matrix element 
is estimated to be of order 0.45 eV leading to a bare 
bandwidth of 3.6 eV. While in principle one expects that 
the Coulomb correlations act to reduce the bandwidth, it 
is not clear what the precise value of the effective band- 
width should be in an RPA framework. Hence, it is nec- 
essary to check the dependence of the results on the ef- 
fective bandwidth W ^ which was assumed to be 1 eV in 
Section HI. For this reason, here, the Knight shifts are 
calculated for = 3 eV. In this case, larger values of U 
are required for fitting the T^^ data and, consequently, 
the pure system has a stronger Stoner enhancement of 
the antiferromagnetic correlations. For W — Z eV, the 
value of the *^Y hyperfine coupling D was chosen to be 
—2 kOe/^B in order to fit the value of ^^K for pure 
YBa2Cu307 in the normal state. The experimental esti- 
mate for D is also about —2 kOe/^s [B. 



In Fig. 13, versus T is shown for optimally doped 
YBa2Cu306+x calculated using W = 3 eV. Here, the ''Li 
hyperfine couphng C was taken to be 1.4 kOe/fiB- The 
open circles represent ^Kq deduced from the ^^K data 
p8[ on pure YBa2Cu307. The dotted curve represents 
the theoretical results for "^Kq. Here, it is seen that a 
smaller value of Vi is required for fitting the data as 
compared to that for = 1 eV. This is because in this 
case the antiferromagnetic correlations are stronger. 

It is also necessary to compare the value of C used 
here with the experimental estimate. In Ref. C was 
estimated to be 2.4 kOe//ZB. This value was deduced by 
comparing the enhancement of the uniform susceptibil- 
ity Ax with the enhancement of '"if for the underdoped 
YBa2Cu306+a; by assuming that the contribution to Ax 
is arising from the changes in the magnetic correlations 
at only the four nearest-neighbour Cu(2) sites of the im- 
purity. Note, however, that Ax is given by 



Ax X! ^("^j) ^ Xpure(q 0), 



(42) 



where the sum over j is carried over the whole lattice, and 
the Knight shift data and the numerical results presented 
in this paper indicate that the changes induced in k(rj) 
by the impurity is extended and not restricted to the four 
nearest- neighbor sites of the impurity. In fact, by limiting 
the changes in k{rj) to be only at the nearest-neighbor 
sites, one would overestimate Ax. Similarly, when one 
calculates C by using the experimental data on Ax and 
the '^K data, one would overestimate C. For instance, in 
Fig. 11(a), it is seen that at 100 K the enhancement of 
A;(1,0) by the impurity, Afc(l,0), is about 20 states/eV 
for underdoped YBa2Cu306+2;. If Ax is estimated by 
just using the change in fc(l,0), then one would obtain 
Ax = 0.4 states/eV in the dilute limit for an impurity 
concentration of 0.5%. On the other hand, the calculated 
value of Ax at T = 100 K and Vi = -0.15t for under- 
doped YBa2Cu306+a: in this model is ~ 0.2 states/eV. 
Hence, by using AA;(1,0) and assuming that the changes 
in fc(r) occur only at the nearest-neighbor sites of the im- 
purity, one would overestimate Ax by about a factor of 
two. Similarly, if AA:(1,0) and Ax are used to calculate 
C one would overestimate C by the same amount. When 
this is taken into account, it is seen that the values of C 
used here, 0.85 kOe/^B in Section HI and 1.4 kOe/fiB 
in Section IV. C, are comparable to the experimental es- 
timate which was obtained by using the data on Ax and 

In Figs. 14(a)-(c), the T dependence of 
—S^^K {outer), and — (5*^i4r(middle) are shown for under- 
doped YBa2Cu306-(-a;. Here, it is seen that for Vi = 
— 0.03t, the magnitudes of the calculated Knight shifts 
agree with the experimental data. However, the fitting 
of the T dependence of the Knight shifts is not as good 
as that seen in Section HI for W = 1 eV. 
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The results presented in this section show that the 
value of Vi required for fitting the Knight shift data de- 
pends on the effective bandwidth. For W = 3 eV, the 
system has stronger enhancement of the AF correlations 
compared to that for W = 1 eV. Here, it has been shown 
that as W increases Vi required for fitting the data de- 
creases. However, even if a W of 4 eV is used, at the 
level of RPA, an impurity potential with Vi = does not 
produce sufficient T variation for fitting the data on 
optimally doped YBa2Cu306+a;, and a weakly attractive 
Vi is necessary. 

D. Nature of the effective impurity potential 

Because of the Coulomb correlations, the bare im- 
purity potential which has only an onsite component 
Vq cJg.C0(T could acquire an extended component 
through higher-order scattering processes ||2^. How this 
could happen at lowest order in the Coulomb repulsion 
U is discussed in Appendix 11. This extended component 
is modelled here by using a static Vi as a free parameter. 
In this picture, the fact that the Knight shift data imply 
14 < is giving information on the effective particle- 
particle interaction in the system. The real-space struc- 
ture of the effective impurity interaction was studied at 
half- filling for the t-J model within the exact diagonal- 
ization calculations It is useful to calculate Vi in 

exact numerical calculations away from half-filling. This 
would be a test of the model presented here. Further- 
more, note that here Vi is assumed to have no temper- 
ature dependence. However, Vi could depend on T, if 
indeed the Coulomb correlations play a role in inducing 
the extended component of the impurity potential. This 
could quantitatively affect the fits seen above. 

These results carried out at the level of RPA imply 
that Vi is weakly attractive. At this level, even though 
an onsite scattering potential (Vi = 0) yields large en- 
hancements of the Knight shifts, it is still insufficient for 
fitting the T dependence, for instance, of in opti- 
mally doped YBa2Cu306+2:. However, note that correc- 
tions beyond RPA might change this result. For instance, 
the spin-fluctuation self-energy corrections to the single- 
particle Green's functions, when taken into account self- 
consistently, could play an important role. 

Within this model, the Knight shift experiments on 
the cuprates with nonmagnetic impurities are interesting 
especially because they probe the interplay of the corre- 
lations in the density and the magnetic channels. The 
effective impurity potential acts in the density channel 
while its spatially-resolved magnetic response is detected 
through the Knight shift measurements. Similarly, the 
recent ^Li Knight shift measurements ]4^ ] in the d-wave 
superconducting state of YBa2Cu306+a; give valuable in- 
formation about the interplay of the density, magnetic 



and the d-wave superconducting correlations in this ma- 
terial. 

Another set of experiments which produce spatially re- 
solved information on the effects of the nonmagnetic im- 
purities are the STM measurements on Zn substituted 
BISCO Q]. These measurements have been carried out 
in the superconducting state. Clearly, similar STM mea- 
surements above Tc would be useful for understanding 
the effects of the nonmagnetic impurities in the normal 
state. This would allow for a direct comparison of the 
single-particle properties with this model. 

It is also desirable to extend these model calculations 
to the superconducting state in order to make com- 
parisons with the STM [||| and the recent NMR 0] 
measurements below Tc. However, note that if a spin- 
fluctuation mediated mechanism is assumed for the d- 
wave pairing, then the effects of the impurity on the 
pairing potential need to be taken into account as well 
as the scattering of the quasiparticles by the impurity. 
This is because it is already known in the normal state 
that in the local environment of the impurity the spin- 
fluctuations are strongly modifled jlj-jl^ . 

Here, the effects of the scatterings by a nonmagnetic 
impurity on the magnetic spectrum is studied by mod- 
elling the impurity as a potential scatterer. In this re- 
spect, disorder could have effects similar to those of non- 
magnetic impurities, and it has been already pointed 
out that the anomalous line broadening of ®^Cu(2) in 
La2-a;Sra;Cu04 could be due to the intrinsic disorder in 
this compound 

V. SUMMARY AND CONCLUSIONS 

In this paper, the Knight shift data on YBa2Cu306+a; 
with nonmagnetic impurities have been analyzed within 
a rather simple model exhibiting short-range antiferro- 
magnetic correlations. The antiferromagnetic correla- 
tions have been modelled within the framework of the 
2D Hubbard model, and the effects of an impurity have 
been approximated by using a static extended impurity 
potential. The strength of the antiferromagnetic correla- 
tions in the pure system has been determined by fitting 
the data of pure YBa2Cu306+a;- The impurity po- 
tential has been assumed to have a range of just one lat- 
tice spacing without any temperature dependence. The 
onsite component of the impurity potential was taken 
to be strongly attractive, and the near-neighbor compo- 
nent Vi was treated as a free parameter in fitting the 
Knight shift data. The simplicity of the model and the 
differences in the stoichiometry of the samples along with 
the uncertainities in the hyperfine couplings are factors 
which limit the conclusions which can be drawn. Nev- 
ertheless, here, it has been found that the anomalous 
T dependence of in YBa2Cu307 can be fitted, if a 
weakly attractive Vi is used. The nature of the effective 
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impurity potential could be studied in the paramagnetic 
state of the Hubbard model, where there are short range 
antiferromagnetic correlations, by using exact numerical 
methods. This could be a test of one of the main as- 
sumptions of this model. Especially, the sign of Vi could 
be tested. These calculations have been also extended to 
the case of underdoped YBa2Cu306+2: after making var- 
ious assumptions about the magnetic correlations in this 
compound, which has a magnetic pseudogap. Based on 
these assumptions, it has been found that the real-space 
structure of the magnetic correlations in the vicinity of 
the impurity is consistent with the Knight shift experi- 
ments. However, caution is necessary in interpreting the 
results on underdoped YBa2Cu306+£!;. 

The results presented in this paper depend on the 
nature of the effective impurity potential used. Any 
magnetic scattering component in VcS could significantly 
change the results. Furthermore, it needs to be kept in 
mind that the validity of these results depend on the 
weak-coupling approach used for calculating the mag- 
netic correlations. Clearly, much remains to be under- 
stood about the effects of the Zn or Li impurities on 
the magnetic correlations of YBa2Cu306+2; in the nor- 
mal state. 



Appendix I: ^Li Knight shift 

In the presence of a uniform static magnetic field, 
B = Bqz, and a hyperfine coupling to the electronic spins 
given by Eq. (nS) , the interaction energy of a ""Li nuclear 
moment fi — h 7„I is 



(43) 



1=1 



where i sums over the four Cu(2) sites neighboring the 
^Li impurity. This expression can be rewritten as 

-h^j„I'Boil+^ K), (44) 

where the ^Li Knight shift ^ii' defined by 

4 



C 



(45) 



gives the fractional change in the Zeeman frequency of 
the nuclear magnetic moment due to the hyperfine cou- 
pling. According to the Kubo linear-response theory Q, 
the expectation value at time t of the electronic spin at 
site Vi, {S''{ri,t)), is given by 



{S'{r,,t))^^i / dt' {[S'-{r,,t),Vz{t')])0{t-t'), 



where Vz(t') = e*^°*Vze with the Zeeman term for 
the electronic spins 



Vz 



(47) 



In Eq. ([461), the expectation values are evaluated with 
respect to the eigenstates of 

Hq = H Hubbard + Vimp, (48) 

where Hnubbard is given by Eq. (£[) and the impurity 
interaction Vimp is given by Eq. (M). In the adiabatic 
limit and as f — > oo, {S^{ri)) is related to the magnetic 
susceptibility Bs], and in this limit Eq. (Efl) reduces to 



{S^{v,))^h^,BoJ2xir^,r,) 



(49) 



where the transverse sus- 

ceptibility x(rj,rj) = x{ri,rj,iuj,n = 0) is defined by 
Eq. (p^). Hence, the ^Li Knight shift is given by 



'In 



where 



(50) 



(51) 



with j summing over the whole lattice. Note that for the 
pure system k(ri) reduces to Xpurei^ 0). 



Appendix II: Ven at lowest order in U 

It is possible that the effective impurity potential has 
an extended component because of the correlated nature 
of the host material ||2^ . At lowest order in the bare im- 
purity potential Vg and the Coluomb repulsion U, such a 
contribution originates from the scattering process shown 
in Fig. 15. This process leads to a momentum-dependent 
effective interaction which is given by 

V}^\q) ^ -VoUx^cO, (52) 

where Xo is the Lindhard susceptibility of the pure sys- 
tem, Eq. (H). Since Xoil) peaks at q ~ (tt, tt), v}^^ is 
attractive at sites neighboring the impurity when Vq < 0. 
For large Vq, it is necessary to replace the bare impurity 
potential in Fig. 15 by the impurity scattering i-matrix, 
in which case one obtains 



(46) 



l-VoFoiiUJn'Y 



(53) 
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where 



Foiiujn) = — 2_^Go(p,iw„). 



The Fourier transform 



(54) 



(55) 



gives the real-space structure of the effective interaction. 
This calculation of T^g^' (r) has been carried out for Vq = 
-lOOi, T = 100 K and = 3 eV, and it is found that 
V}g\r) at r = (1,0) is -0.017i. Even though v}g\r = 
(1,0)) is found to be attractive in leading order in U, 
its magnitude is smaller by a factor of three from Vi = 
— 0.05t found by fitting the "^K data on optimally doped 
YBa2Cu306+a;- 

This is the behavior expected for Vcs at lowest order 
in U. Clearly, an approximation which is first order in U 
would be insufficient. The actual structure of 14^, which 
includes scattering processes to all orders in U, depends 
on the reducible particle-particle vertex in the singlet and 
the triplet channels. Hence, it would be useful to calcu- 
late VcS using exact numerical techniques in the param- 
agnetic state of the Hubbard or the t-J models. However, 
it should be noted that the structure in the electron den- 
sity n(r) does not necessarily reflect the structure in VoS, 
as it was seen in Fig. 11(c). In Ref. [^, Kff was calcu- 
lated for the t-J model but at half-filling, where VcS at 
the (1,0) site was found to be repulsive. However, the 
insulating state might be different than the paramagnetic 
state, which is considered here. In Ref. an effective 
impurity interaction which is repulsive at r = (1,0) was 
used for calculating the enhancement of the uniform sus- 
ceptibility by dilute nonmagnetic impurities. Since both 
positive and negative values of Vi act to enhance the 
uniform susceptibility, this does not determine the sign 
of Vi. On the other hand, the Knight shift data studied 
here require that, at the level of RPA, VeS at the nearest- 
neighbor sites of the impurity is weakly attractive. 
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FIG. 1. Illustration of the extended impurity poten- 
tial due to a nonmagnetic impurity represented by the filled / . \ 

square at the center. Here, Vb is the onsite component of the ^ ' 

impurity potential and Vi acts at the nearest-neighbour sites 

of the impurity. 1.5 ^ '''' '''''^ '''''''' ' 

100 200 300 400 

T (K) 

FIG. 2. (a) Transverse nuclear relaxation rate T^^ for 
®^Cu(2) versus T in pure YBagCusOe+a;. The open and 
the filled circles represent the experimental data by Imai 
et al. [24] and Takigawa [25] for YBa2Cu306.9 and for 
YBa2Cu3 06.63, respectively. The dashed line is a fit of 
the Tj"^ data by Imai et al. [24] on YBa2Cu306.9. For 
YBa2Cu306.63, the fitting of the data below 160 K has been 
carried out in two ways. The solid line has been obtained 
by ignoring the saturation of T^^ below ~ 160 K, and 
the dotted curve has been obtained by fitting the data for 
80 K< T < 160 K. (b) Temperature dependence of U. The 
dashed curve has been obtained by fitting the T^^ measure- 
ments on pure YBa2Cu30(i.9. The solid and the dotted curves 
have been obtained for pure YBa2Cu3 06.63 by fitting the solid 
and the dotted curves seen in (a). 
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FIG. 3. (a) XpuT-e(q) versus q at 100 K for the pure sys- 
tem. The dashed curve represents results for optimally doped 

YBa2Cu306+a;. The solid and the dotted curves have been ob- 
tained for underdopod YBa2Cu306+x. (b) x(Q*) versus T for 
optimally doped YBa2Cu306+x, where Q* is the wavevector 
at which x(q) peaks at low temperatures. Here, the dashed 
curve is for optimally doped pure YBa2Cu306+a;, and the 
solid curve is for the case of 0.5% dilute impurities, (c) x(Q*) 
versus T for underdoped YBa2Cu306+a:. The dashed curve 
is for the pure system and the solid curve is for the case with 
0.5% impurities. The dotted curves were obtained by using 
the values of U given by the dotted curves in Fig. 2(b). The 
results for the impure systems shown in (b) and (c) were cal- 
culated for Vi = — 0.15t as described in Section IV. A. 
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FIG. 5. Sketch of the neighbouring Y(i) sites of a non- 
magnetic impurity located in the Cu02 plane. Here, the filled 
square indicates the impurity site. 



FIG. 4. Feynman diagrams for (a) the dressed sin- 
gle-particle Green's function G and (b)-(c) the irreducible 
susceptibility xo within the presence of one impurity. 



18 




E 

Q. 
Q. 



100 - 



I I I 



I I I I 



(b) 



100 



200 
T (K) 



300 



400 



FIG. 6. (a) ^Li Knight shift "^K versus T in optimally 
doped YBa2Cu306+a:- The filled circles represent the Knight 
shift data by Bobroff et al. [2]. The curves represent the re- 
sults of the calculations for various values of Vi/t which are 
indicated next to the curves, (b) Estimated temperature de- 
pendence of ^ Kq for optimally doped YBa2Cu306-i-a!. Here, 
Kq represents the value of the ^Li Knight shift, if the sub- 
stitution of the Li impurity had not induced any changes in 
the magnetic correlations around it. The dashed curve is the 
experimental estimate of ''Kq obtained by using the ^^K data 
[38] on pure optimally doped YBa2Cu306+x- The solid line 
is the result obtained by setting Vb = Vi = in this model. 
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FIG. 7. (a) Temperature dependence of —5 K at the sites 
which are the 1st, 2nd, 3rd and the 4th nearest ^^Y neigh- 
bors of the impurity in optimally doped YBa2Cu3 06+i^. Here, 
5^^'' K is the induced change in the *^'*Y Knight shift up on the 
substitution of the impurity, (b) Temperature dependence of 
the *'^Cu(2) Knight shift for H||c, "^ATcCri), at various sites 
Vi, which are indicated next to the curves, for the optimally 
doped YBa2Cu306+a:. The results shown in (a) and (b) were 
obtained by using Vi = — 0.15t. 
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FIG. 8. (a) Temperature dependence of 5^ K for under- 
doped YBa2Cu306+2!. Here, 5'' K is defined as ^ K —"^ Ko, 
whore is the ^Li Knight shift and '^Ko represents the 
value of the ^Li Knight shift, if the substitution of the Li 
impurity had not changed the magnetic correlations around 
it. The filled circles represent the data by Bobroff et al. [2] on 
YBa2Cu306.6. The curves represent the results of the calcula- 
tions for various values of Vi /t which are indicated next to the 
curves, (b) Temperature dependence of ^Ko for underdoped 
YBa2Cu306+x. The solid points have been obtained by us- 
ing the K{ma,m) data on underdoped YBa2Cu306+a; with 
Zn impurities [1,4], while the dashed curve has been obtained 
by using the '^'^ K data on pure underdoped YBa2Cu:iO(i+,r 
[38] . The solid curve shows the results of the calculations for 
Vo = Vi= 0. 
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FIG. 9. (a) Temperature dependence of — (outer) 
in underdoped YBa2Cu306+a! with Zn impurities. Here, 
S^^ K{outci) is defined as the shift of the outer '^"Y satellite 
with respect to the ^''Y main resonance line. The open and the 
filled circles represent the experimental data on — 5*®_R'(outer) 
in underdoped YBa2Cu306+a; with Zn impurities for the mag- 
netic field H]]c and H ± c, respectively [1,4]. The curves 
represent the results obtained by using various values of Vi/t 
which are indicated next to the curves for the Y(l) site, which 
is the nearest ^®Y neighbour of the nonmagnetic impurity, (b) 
Results similar to those in (a) but for — 5*® (middle). Here, 
5^^ K (middle) is defined as the shift of the middle ®®Y reso- 
nance line with respect to the main line. The curves represent 
results for various values oiVi/t for the Y(2) site, which is the 
second-nearest ^^Y neighbour of the nonmagnetic impurity. 
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FIG. 10. Temperature dependence of the ''^Cu(2) Knight 
shift for H||c, ^^Kc{ri), in underdoped YBa2Cu306+2! at sites 
Vi which arc indicated next to the curves. These results were 
obtained for Vi = — O.lSt using the values of U{T) given by 
the solid curve in Fig. 2(b). 
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FIG. 11. (a) fc(r) versus the distance r away from the 
impurity in units of the lattice spacing o for T = 100 K and 
V\ — — 0.15t. Here, the open and the filled circles repre- 
sent the results for the optimally doped and the underdoped 
YBa2Cu30(i+a;, respectively. Also, the horizontal long-dashed 
and dotted lines represent k{r) for the optimally doped and 
underdoped pure systems, respectively, (b) x(r, r) versus the 
distance r away from the impurity presented in the same way 
as in (a) for fe(r). (c) Electron occupation n{ri) at sites rj 
near the impurity plotted as a function r = \ri\ at 100 K. 
These results were obtained for Vi/t = —0.15 with W = 1 eV 
(solid) and for Vi/t = -0.05 with = 3 eV (dashed). Also 
shown are results for an onsite impurity potential (Vi = 0) 
and W = 1 eV (dotted). In all of these cases, Vb = — lOOt 
was used. 
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FIG. 13. Temperature dependence of for optimally 
doped YBa2Cu306+a;, obtained for various values of Vi/t. 
These results arc for an effective bandwidth of 3 oV, and here 
the ''Li hyperfine coupling C was taken to be 1.4 kOe/fiB- 
The filled circles represent the ^ K data on optimally doped 
YBa2Cu306+x by Bobroff et al. [2]. The open circles rep- 
resent the estimate of Kg obtained by using the ^^K data 
[38] on pure optimally doped YBa2Cu306+a^. Here, the 
hyperfine coupling was taken to be -2 kOe//is. The dotted 
curve has been obtained by using Vb = Vi = and the re- 
maining curves have been obtained by using Vb = — WOt and 
the values of Vi jt which are indicated next to the curves. 



FIG. 12. Off-diagonal susceptibilities (a) — Xo(q, q' = 0) 
and (b) — xCQ: q' = 0) versus q obtained by using Vx = — 0.15t 

for optimally doped YBa2Cu306+:r at T = 100 K. Here, the 
(5-function components at q = q' = have been omitted. 
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FIG. 14. Temperature dependence of (a) S^K, 

(b) -5^^K{outer) and (c) -5*'^ 7f (middle) for underdoped 
YBa2Cu306+a; obtained by using various values of Vi/t which 
are indicated next to the curves. These results have been 
calculated by using an effective bandwidth W oi 3 eV, as 
compared to W = 1 eV used in Figs. 8 and 9. In addition, 
here the '^Li and the hyperfine couplings were taken to be 
1.4 kOo//is and -2 kOo//is, respectively. The experimental 
data in (a) are from Ref. [2], and the data in (b) and (c) are 
from Refs. [1,4]. 
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FIG. 15. Feynman diagram for a scattering process which 
causes the effective impurity potential to be extended in real 
space. This diagram is lowest order in U and the bare im- 
purity potential Vo. Here, an electron with momentum p, 
Matsubara frequency w„ and spin a scatters to a state with 
momentum p + q, Matsubara frequency w„ and spin a. 
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